Kinematical Test of Large Extra Dimension in Beta Decay Experiments 
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The forthcoming experiments on neutrino mass measurement using beta decay, open a new window 
to explore the Large Extra Dimension model. The Kaluza-Klein tower of neutrinos in Large Extra 
Dimension contribute to the Kurie function of beta decay that can be tested kinematically. In 
addition to providing an alternative approach using just the kinematical properties, we show that 
KATRIN can probe the compactification radius of extra dimensions down to 0.2 /im which is better, 
at least by a factor of two, than the upper limits from neutrino oscillation experiments. 
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The "hierarchy problem" is one of the long stand- 
ing problems of the Standard Model (SM) of particles 
which triggered several beyond standard model theo- 
ries (such as supersymmetry, technicolor, . . .) as a way 
out. As a simple description, the hierarchy problem is 
the large disparity between the weak interaction scale 
Mew ~ 10 3 GeV and the scale at which the gravity 
is strong, the so-called Planck scale Mpi ~ 10 19 GeV. 
Within the SM particles, the anticipated Higgs particle 
suffers from this hierarchy through the radiative correc- 
tions to its mass which result in a Higgs mass of the 
order of Planck mass, contrary to the expectation from 
electroweak precision tests. 

The theory of Large Extra Dimensions (LED) pro- 
posed in the seminal papers [H, HJ as an elegant solu- 
tion to the hierarchy problem. The basic assumption in 
LED is that there is just one fundamental scale, which 
is the weak scale Mew, an d the large value of the ob- 
served Mp\ from gravity is a manifestation of the exis- 
tence of extra dimensions such that Mpi ~ Mew in the 
higher dimensional space but effectively is large in our 
4-dimcnsional space. Assuming the existence of n com- 
pactified extra dimensions, it can be shown that gravity 
would deviate from the inverse-square law at distances 
~ 10ir- 17 (TcV/M EW ) 1+ " cm Q. Thus, to avoid the 
conflict with the confirmed inverse-square behavior of 
gravity at solar scale, the number of extra dimensions 
should be n > 2. However, here we assume an asymmet- 
ric space such that one out of the n extra dimensions is 
compactified on a larger spatial scale -Red and effectively 
we are facing a 5-dimcnsional space. 

The setup of LED model is as follows: the particles 
with charge under the SM gauge group including the 
charged leptons, active neutrinos, quarks, gauge bosons 
and Higgs particle live on a 4-dimensional brane embed- 
ded in the (4 + n)-dimcnsional space. Several mechanism 
can justify the localization of these particles [H, H|, but to 
justify Mpi <~ Mew (weakness of gravity) it is assumed 
that the singlets of the SM gauge group including the me- 
diator of gravity (graviton) can propagate freely in the 
extra dimensions. 

A bonus of the LED model is to give a natural explana- 
tion for the smallness of neutrino masses 0-0] • As a con- 



crete model, let us consider the LED extension of the SM 
augmented by three masslcss right-handed 5-dimcnsional 
neutrinos \1/ Q corresponding to the three active neutrino 
flavors. To be singlet under the SM gauge group, en- 
able the $ a fields to propagate in the 5-dimensional 
(x 11 ,'}/) space, where y denotes the large extra dimen- 
sion (/x = 0,1,2,3). The general action describing the 
neutrino sector (assuming lepton number conservation by 
assigning lepton number -1 to is given by 

S = j d 4 xdyzf Q r A a A * Q + J d A x {ivS'y^vS 
+ K^HP^(x,y = 0)+-R.c.}, (1) 

where A = 0, . . . , 4; H is the Higgs doublet, k is the 
Yukawa coupling matrix and the right-handed neutrino 
field is decomposed to (V'liV'r)- After electroweak sym- 
metry breaking, the above action results in the 3x3 
Dirac neutrino mass matrix ~ kv / ' y/V n M^jj, where v is 
the VEV of Higgs doublet and V n is the volume of the 
extra n-dimensional space which naturally suppresses the 
mass matrix. Rotating the neutrino states to the mass 
eigenstates, gives the three eigenvalues of the 3x3 mass 
matrix (mf, , m§). 

From the 4-dimcnsional brane point of view, the right- 
handed neutrinos appear as a tower of Kaluza-Klein 
(KK) modes with increasing masses. Taking into account 
the quantized Dirac masse of these KK modes, the eigen- 
values of the resulting infinite-dimensional mass matrix 
gives the masses of the neutrino mass eigenstates. The 
eigenvalues are given by Xi/Rei) (i = 1, 2, 3), where Aj is 
the root of the following transcendental equation 043 



Aj — 7ri?ED ( m ?)~ cot(7rAi) = 



(2) 



It is easy to check that Eq. ([2]) has infinite number of so- 
lutions X- , where n = 0, 1, . . .; and X^ G [n, n+1/2]. It 
should be emphasized that (mf, m!?, m^) are the param- 
eters of the action in Eq. ([T]) and the masses of the neu- 
trino states (including the three active SM neutrinos and 
the KK tower of sterile neutrinos) are given by \[ n ^ / Red 
(i = 1,2,3), where n = and n = 1, . . . , corresponds 
respectively to the active and sterile neutrinos. However, 
the two solar and atmospheric mass-squared differences 
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(Am 2 ol , Am 2 tm ), from neutrino oscillation phenomenol- 
ogy, fixes two out of the three parameters (m® , m§ 7 m^). 
The procedure is as follows: regarding and -Red as 
free parameters, it is possible calculate the value of A^ 
from Eq. ^ for fixed values of mf and -Red- Then, the 



values of A 2 and X^' is given by 



(0) 



and 



(a|°)) 2 = (a^) : 



£ ED Am 2 tm , 



(3) 
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for the normal hierarchy (NH) scheme on neutrino masses 
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scheme (A^ 
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a 2 > A x 
free parameters, A 3 



For the inverted hierarchy (IH) 
> A, ), regarding and -Red as 
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^ can be calculated from Eq. ([2|), and 



10 u 



> 10 

Q 
LU 

10 2 



10'- 



10 



10" 



10" 



10"' 



I I 1 1 

" R ED = 10" 7 m 

- 

r 




- 

- 
: 






M"' ' ^ED 






' ^ED 




s 


4 0) /R ED 





10" 



10 u 



10' 



mrteV] 



FIG. 1. The plot of \f ) /R ED (i = 1, 2, 3) as a function of mf 
assuming NH scheme and Red = 10 -7 m. 
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Having the values of A^ and Ag "* for NH (A^ and A^ 
for IH), it is possible to read the value of and 
(or mf and for IH) from the Eq. ([2]). After that, the 
masses of sterile states in the KK towers can be calculated 
by finding the roots A,-™^ "* of Eq. ([2]). Thus, the only free 
parameters in the mass matrix of the action Eq. (Q~J) are 
the (m?,i?ED) for NH and {m$,R ED ) for IH. 



However, in deriving the values of A 2 °^ and A 3 u; from 
Eqs. ()3I4|) . or A^ and A 2 °^ from Eq. (0 in the case of 
IH, we are constraint by the fact that for any values 
of (Am 2 ol , Am 2 ^) always we should have x[ ^ < 0.5 
(see the comment after Eq. ([2])). This inequality de- 
fines the allowed physical region in the (hi^-Red) or 
(m^, -Red) parameter space which a solution exist to the 
observed solar and atmospheric oscillation scales within 
the framework of LED model. Fig. [TJ shows the behavior 
of A,|°^/-Red as a function of assuming NH scheme 
for (Am s 2 ol , Am 2 tm ) = (7.6 x 1(T 5 ,2.4 x 10" 3 ) eV 2 [| 
and -Red = 10~ 7 m. As can be seen, for <C 1/-Red 
we have A^ /-Red — which comes from the fact that 
in this case the KK tower sterile neutrinos have a very 
small contribution to the neutrino mass matrix. But, 
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however, for 



~ 1/-Red the mass A^V-Red saturate 
to the value 0.5 /-Red ■ The black vertical thick line in 
Fig. [T] shows the discussed allowed region coming from 
the condition A^ < 0.5 . 

The contribution of the KK tower sterile neutrinos to 
the neutrino mass matrix leads to the mixing between the 
active and sterile neutrinos. The active flavor neutrino 
states can be written as 0-0] 



i=l 



n=0 



(6) 



where are the mass eigenstates composed of the n th 
KK mode of iph (n — corresponds to the SM active mass 
eigenstates) and the 3x3 matrix U is the Pontecorvo- 
Maki-Nakagawa-Sakata mixing matrix of neutrinos 0]. 
The parameters L\ n come from the diagonalization of 
the infinite-dimensional mass matrix |5j-l7j and is given 
by 



{L 0n f = 



1 + tt2 (R EB mf) +XY l)2 /(R ED 7 



(7) 



The values of the L\ n parameters also is just a func- 
tion of the two free parameters (to^-Red) for NH 
and (t^j-Red) for IH. To clarify the behavior of the 
parameters in Eq. ([7]) we plotted in Fig. [2] the L?" 
for n — 0, . . . , 9; as a function of and assuming 
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FIG. 2. The plot of L? n as a function of m" assuming 
Red = l(T 7 m. The curves from top to bottom correspond 
to n = 0, . . . , 9. 
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the values of L® 71 ^ is very small which means that the 
contribution of the sterile mass eigenstates to the decom- 
position of active flavor states in Eq. ^ is almost negligi- 
ble. But, for mf ~ l/i?ED, KK tower states with higher 
masses play role in the active states decomposition. 

The decomposition of the flavor neutrino states in 
Eq. ([6]) would impact the phenomenology of the neutrino 
sector of SM in a substantial way. From flavor oscillation 
point of view, considering the decomposition in Eq. ([6]), 
the probability of — > oscillation for neutrinos with 
energy E v at the baseline L is given by 0, 0] 



P, 



a/3 



3 oo 

J2u ai U^*Y^{Li n ) 2 exp 

i = l n=0 
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(8) 



Plenty of works is devoted to use of the above for- 
mula in fitting the data of oscillation experiments @, 
an d recently 12] . The updated result includ- 
ing CHOOZ, MINOS and KamLAND experiments gives 
" 7 m at 90% C. 



12] 



Red < 6 x 1CT 7 m at 90% C. L. 

From the kinematical point of view, the decomposi- 
tion in Eq. ([6]) would impact the neutrino mass mea- 
surement experiments. Specifically, the KK tower sterile 
neutrino states contribute to the effective electron anti- 
neutrino mass which can be probed at the beta decay 
experiments. The forthcoming KArlsruhe TRItium Neu- 
trino (KATRIN) experiment |l3(, probes the kinemat- 
ical effect of the D e produced in the tritium beta de- 
cay 3 H— > 3 Hc + + e~ + D e by the measurement of elec- 
tron energy spectrum near the endpoint of the reaction 
Q = 18571.8 ± 1.2 eV Q. In the absence of neutrino 
mass, the Kuric function of beta decay is a linear function 
of the electron's kinetic energy T e , that is K(T e ) = Q—T e . 
A nonzero value of the neutrino mass leads to distortion 
in the Kurie function, such that the endpoint of the elec- 
tron's energy spectrum displaces by the lightest neutrino 
mass value and heavier masses generate kinks in the spec- 
trum. Assuming the decomposition of Eq. (|6]) for electron 
anti-neutrino, the Kurie function takes the following form 

K(T e ,m ,R E v) = (9) 



n=0 
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°k 



X 



(n) 



-Red 



where £k = Q — Wk — T e ; and hereafter mo = mf for NH 
and mo = m^ for IH. The dependence on mo comes from 
the \\ nS> (see Eq. ([2])). The Wk and pk are respectively 
the excitation energy and transition probability of the k th 
rotational and vibrational excited state of the daughter 
nucleus in the tritium beta decay (loT ]. Also, inserting 
the Heaviside step function 0(£k — A] /Red) in Eq. © 
guarantees the conservation of energy. 

The spectrometer of the KATRIN experiment is accu- 
mulative, which means that it can collect the electrons 
with kinetic energy larger than the electrostatic barrier 



qU which can be tuned. The rate of the electrons passing 
the potential barrier is given by 

S(qU,m ,R EB ) = N b + (10) 

/ F(Z, T e )E ePe K(T e , m , i? ED )i?'(T e , qU) dT e , 
Jo 

where E e and p e are respectively the energy and momen- 
tum of the electron. The F(Z, T e ) is the Fermi function 
which takes into account the electromagnetic interaction 
of the emitted electron in beta decay with the daughter 
nucleus (Z = 2) [l6j]. The function R' shows the trans- 
mission function of the KATRIN experiment including 
the resolution of the spectrometer and the energy loss 
processes in the propagation of electrons from the source 
to collector. Finally, the A^ s is the normalization factor 
calibrating the rate of electron emission which for KA- 
TRIN is 1.47 x 10~ 13 s" 1 eV" 5 ; and N b is the background 
rate which is 10 mHz jl3j|. The details of the KATRIN 
experiment and functional form of R' can found in [l7| . 

We calculated the sensitivity of KATRIN experiment, 
after three years of data-taking, using the proposed op- 
timized running time for a barrier potential qU £ [Q — 
20,Q + 5] eV (see (H). Fig. |3] shows the 90% C. L. sensi- 
tivity to the (mo, -Red) parameters assuming NH for the 
neutrino mass scheme. The black solid line shows the 
boundary of physical region of parameters which comes 
from the requirement that A| ' < 0.5 ; i.e., the values 
of (mo, -Red) parameters in the right-hand side of the 
black curve is not compatible with (Am s 2 ol ,Amj tm ) = 
(7.6 x 10~ 5 , 2.4 x 10~ 3 ) eV 2 . The red solid curve shows 
the 90% C. L. sensitivity of KATRIN (for NH). To com- 
pare with the limits from oscillation experiments, we cal- 
culated the limit from CHOOZ experiment (blue dashed- 
dotted curve), which is the strongest limit within oscil- 
lation experiments. For the CHOOZ experiment [l8| . 
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FIG. 3. The 90% C. L. sensitivity in (mo, Red) parameter 
space for KATRIN experiment assuming NH for the neutrino 
mass scheme. Also we show the 90% C. L. limit form CHOOZ 
experiment and the combined analysis of KATRIN+CHOOZ. 
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FIG. 4. The same as Fig. [3] but for IH for neutrino mass 
scheme. Also, in this case the combined KATRIN+CHOOZ 
is the same as KATRIN alone. 



which is a reactor v e disappearance experiment at the 
base-line ~ 1 km, we fit the obtained data using Eq. ((HJ 
for a = f3 = e. As can be seen, the KATRIN sensitivity is 
about one order of magnitude stronger in -Red ^ 10~ 7 m 
region; but weaker for -Red ^ 10 -6 m. To take advantage 
of the oscillation experiments, we did the combined anal- 
ysis of KATRIN+CHOOZ which is shown by the green 
dashed line. 

In Fig. Hwe show the 90% C. L. sensitivity of KATRIN 
for the IH case. As can be seen, for all the values of too 
and -Red, the KATRIN sensitivity is stronger than the 
CHOOZ limit which is the strongest within the oscillation 
experiments. Also, because of the stronger sensitivity in 
KATRIN, the combined analysis of CHOOZ+KATRIN 
gives the same result as KATRIN alone. 

The KATRIN sensitivity is much stronger than the 
limits coming from probing deviations of gravity law 
from inverse-square at sub-millimeter scales, which is 
-Red < 1-6 x 10~ 4 m [l9[. On the other hand, strong 
limits on the size of extra dimensions exist from the as- 
trophysical and cosmological considerations. The astro- 
physical limits stem from cooling time of supernovae [20j | 
and diffuse gamma rays from the decay of KK gravi- 
tons produced at the supernovae For the number 

of extra dimensions n = 2 (n = 3), the upper limit 
from cooling of supernovae is -Red < 0.66 x 10~ 6 m 
(-Red < 0.8 x 10~ 9 m) and the limit from KK gravi- 
tons decay to photons is -Red < 3.61 x 10~ 7 m (-Red < 
3.95 x 10~ 10 m). The production of KK modes in the 
early universe also can significantly affect the cosmologi- 
cal evolution (22[ ■ The limits from this cosmological con- 
sideration is i?ED < 2.2 x 10~ 8 m (-Red < 2.5 x 10" 11 m) 
for n = 2 {n = 3). However, it should be noticed that 
all of these limits from astrphysics and cosmology suffer 
from model dependency. For example, it can be shown 
that compactification of extra dimensions on a hyperbolic 



manifold relaxes all the above mentioned limits com- 
pletely Q . Specifically, for the model we are considering 
in this paper, with an asymmetrical compactification of 
extra dimensions, the limits from astrophysics and cos- 
mology do not apply (see also [24J) The robustness of the 
limit from forthcoming KATRIN experiment is based on 
the fact that limit comes just from kinematical consid- 
erations and do not suffer from theoretical uncertainties. 
In this sense, the KATRIN experiment triggers an al- 
ternative approach to the test of large extra dimension 
model. 

In summary, we have shown that the forthcoming neu- 
trino mass measurement experiments, using the beta de- 
cay (such as KATRIN) , can test the LED model extended 
with singlet 5-dimcnsional neutrino fields that can prop- 
agate in the bulk. This extension was motivated as an 
elegant way to justify the smallness of neutrino masses. 
Besides the fact that this kinematical test at KATRIN 
is a new window to probe LED, we have shown that 
this test can give stronger limit on the compactification 
scale -Red- For the case of inverted hierarchy in the neu- 
trino mass scheme, KATRIN can obtain the upper limit 
-Red < 2.3 x 10" 7 m (90% C. L.) after three years of data- 
taking (for mo — > 0); which is a factor of two better than 
the limit from oscillation experiments. Also, for the nor- 
mal hierarchy scheme, KATRIN can exclude the regions 
of the (too,.Red) parameter space which is inaccessible 
to oscillation experiments. 
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